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I.  EXECUTIVE  SUMMARY 

APL-UW  made  acoustic  measurements  of  near-surface  bubbles  in  fetch-limited 
shallow  water  in  January  and  February  1986  as  part  of  the  High  Frequency  (10  to 
500  kHz)  Acoustics  Project  managed  by  NOARL.  This  work  is  applicable  to  torpedo 
guidance,  control,  and  countermeasure  technologies,  and  particularly  to  evaluating  detec¬ 
tion  and  signal-processing  constraints  imposed  by  a  near-surface  environment. 

Bubbles  near  the  surface  were  ob'  erved  indirectly  by  measuring  vertical  incidence 
backscatter  from  the  near-surface  volume  and  producing  vertical  scattering  strength 
profiles.  Bubbles  were  observable  from  these  measurements  at  wind  speeds  as  low  as 
3  m/s.  Surface  forward  scattering  measurements  at  low  grazing  angles  were  also  taken. 
The  computed  decibel  loss  from  the  forward  scattering  measurements,  attributable  to 
bubbles  and  referred  to  as  surface  bubble  loss  (SBL),  was  greater  than  10  dB  for  wind 
speeds  greater  than  about  5  m/s.  Direct  estimates  of  the  integrated  volume  scattering 
strength  O/  were  obtained  from  vertical  incidence  measurements  (VI  data),  and  inferred 
estimates  of  Oj  were  obtained  from  the  SBL  estimates  derived  from  the  forward  scatter¬ 
ing  loss  measurements  (FL  data).  The  principal  results  are 

(1)  The  depth  dependence  of  the  near-surface  volume  scattering  strength  fol¬ 
lowed  the  expected  exponential  behavior  in  the  majority  of  cases,  but 
significant  departures  occurred  in  some  cases  where  a  two-part  exponential 
behavior  is  observed. 

(2)  The  frequency  dependence  of  O/  derived  from  VI  data  was  o7  «  /2  4;  the  fre¬ 
quency  dependence  of  O/  inferred  from  FL  data  was  Oj  «  /°-6. 

(3)  The  wind  speed  dependence  of  Oj  derived  from  VI  data  was  G/  UbA ;  the 
wind  speed  dependence  of  C/  inferred  from  FL  data  was  o7  «  U22. 

(4)  The  differences  in  the  VI  and  FL  data  noted  in  (2)  and  (3)  may  in  part  be  due 
to  the  VI  data  being  less  sensitive  than  the  FL  data  to  the  near-surface  bubble 
layer. 

The  data  presented  here  represent  a  modest  sampling  of  conditions  in  fetch-limited 
seas.  These  conditions  differ  from  those  in  the  open  ocean,  e.g.,  in  the  length  scales  of 
sea  surface  waves  and  characteristics  of  bubble  plumes,  so  wind  speed  and  frequency 
power  laws  reported  here  may  also  differ  from  those  in  the  open  ocean. 
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II.  INTRODUCTION 

This  report  presents  results  of  acoustic  measurements  of  vertical  incidence  back- 
scattering  and  forward  loss  from  near-surface  bubbles.  The  measurements  were  made  off 
Whidbey  Island  in  Puget  Sound  in  January  and  February  1986  using  1  ms  pulsed  cw  sig¬ 
nals  with  center  frequencies  ranging  from  15-50  kHz.  This  work  is  part  of  the  High- 
Frequency  Acoustics  Project  (administered  by  the  Naval  Oceanographic  and  Atmos¬ 
pheric  Research  Laboratory),  which  primarily  involves  exploratory  development  within 
the  frequency  range  of  10  to  500  kHz.  Developments  are  applied  to  torpedo  guidance 
and  control  and  countermeasure  technologies,  such  as  detection  and  processing  in  a 
near- surface  environment.  A  companion  report  (Dahl  and  McConnell,  1990)  covers  hor¬ 
izontal  and  vertical  spatial  coherence  measurements  made  at  the  same  experimental  site. 
Here,  we  combine  and  extend  earlier  results  from  this  experiment  presented  by  McCon¬ 
nell  (1988a,b). 
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HI.  EXPERIMENTAL  ARRANGEMENT 

Figure  1  shows  the  measurement  location  off  Whidbey  Island  in  Puget  Sound.  We 
chose  this  site  for  its  fairly  open  exposure  to  passing  storms  at  varying  fetches  up  to 
50  km  (for  winds  from  the  west).  Mile  long  cables  connected  the  acoustic  and  environ¬ 
mental  instrumentation  to  a  shore  site  trailer  containing  electronic  systems,  computer, 
and  recorders. 


Figure  1 .  Experiment  site  off  Whidbey  Island  in  Puget  Sound. 


Data  reported  here  were  collected  under  winter  environmental  conditions  character¬ 
ized  by  mean  wind  speeds  ranging  from  0-15  m/s,  with  gusts  recorded  up  to  23  m/s,  and 
rms  wave  heights  ranging  from  2-50  cm.  The  wind  speed  data  originate  from  the 
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National  Ocean  Survey’s  Smith  Island  recording  station  (about  4  miles  from  the  experi¬ 
mental  site)  and  from  a  cup  anemometer  located  on  shore.  The  Smith  Island  data  are  an 
average  over  1  hour  measured  approximately  5  m  from  the  sea  surface,  and  the  cup 
anemometer  data  are  an  average  over  approximately  10  minutes.  All  wind  speeds 
reported  here  have  been  converted  to  the  standard  value  of  10  m  above  the  sea  surface 
using  the  method  of  Smith  (1988).  We  determined  the  rms  surface  wave  height  from  the 
travel  time  fluctuations  produced  by  the  air-sea  interface  return  from  the  vertical 
incidence  backscatter  signal.  Typical  air  temperatures  ranged  from  4°  to  9°C,  with  the 
sea  temperature  near  the  surface  being  a  steady  4.35°C.  The  sound  speed  profile  was 
measured  during  the  course  of  the  experiment  using  CTD  casts.  The  winter  conditions 
caused  strong  mixing  within  the  water  column,  resulting  in  a  nearly  isovelocity  profile  of 
mean  speed  1460  m/s. 

The  measurement  geometry  is  shown  in  Figure  2.  We  made  vertical  incidence 
measurements  using  two  1.3  m  long,  linear  arrays  mounted  perpendicular  to  each  other 
atop  a  5  m  tripod  tower  on  the  seabed.  Acoustic  pulses  were  transmitted  on  one  array 
and  received  on  the  other,  producing  a  very  narrow  combined  beam  with  a  3  dB 
beamwidth  ranging  from  1.2°  at  50  kHz  to  4.1°  at  15  kHz.  The  frame  on  which  the  trans¬ 
ducers  were  mounted  could  be  varied  in  elevation  (0°  to  90°)  and  azimuth  (300°  total). 
The  principal  type  of  measurement  was  vertical  incidence  backscatter  at  15-50  kHz 
using  narrow  beams  and  short  pulse  lengths.  From  these  data  we  determined  the  near¬ 
surface  volume  scattering  strength  profile  attributable  to  bubbles. 

The  volume  scattering  strength  profile  can  be  related  to  the  scattering  strength  for 
an  arbitrary  bistatic  geometry  and  therefore  the  surface  forward  loss.  To  test  this  rela¬ 
tion,  we  interspersed  surface  forward  scattering  loss  measurements  at  low  grazing  angle 
(plus  one  low  angle  backscatter  measurement)  with  the  vertical  incidence  measurements. 
For  the  surface  forward  scattering  measurements  we  placed  a  nearly  omnidirectional 
source  5  m  above  the  bottom  at  a  range  of  500  m  from  the  tower  as  depicted  in  Figure  2. 
Signals  were  received  by  a  nearly  omnidirectional  hydrophone  placed  next  to  the  line 
arrays  on  the  tower.  The  alignment  of  the  acoustic  source  and  receiver  was  roughly 
parallel  to  the  shoreline.  Winds  were  primarily  offshore  or  onshore,  so  the  wave  crests 
were  approximately  parallel  to  the  direction  of  acoustic  propagation.  Because  of  the 
35  m  water  depth,  the  receiver  detected  a  number  of  multipath  arrivals,  including  the  sin¬ 
gle  surface  bounce  path  of  primary  interest.  The  dominant  arrival  involving  the  surface 
is  this  single  bounce  path  since  paths  involving  the  bottom  suffer  a  high  loss  owing  to  the 
porosity  (0.7)  of  the  soft  mud  at  the  measurement  site. 

Both  the  vertical  incidence  backscatter  and  surface  forward  scattering  measure¬ 
ments  used  a  sequence  of  100  short  pulses  for  each  data  set.  For  the  first  part  of  the 
experiment  the  data  were  recorded  digitally  at  a  single  frequency.  Pure  tone  pulses  with 
lengths  of  0.3  ms  and  1  ms  were  interleaved  and  transmitted  at  a  rate  of  one  every  2  s. 
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Figure  2.  Sketch  of  the  measurement  geometry  for  vertical  incidence  backscattering 
measurements  using  two  perpendicular  line  arrays,  and  forward  loss  meas¬ 
urements  using  a  broadbeam  transmitter  and  receiver.  Also  shown  are  the 
physical  processes  responsible  for  producing  and  influencing  near-surface 
buobies  associated  with  wave  breaking. 


giving  a  total  averaging  time  of  200  s  for  50  pings  of  the  same  pulse  length.  The  shorter 
pulse  length  data  were  used  primarily  to  distinguish  multipaths  and  are  not  discussed 
here.  During  the  latter  part  of  the  experiment,  the  digital  data  acquisition  system  became 
inoperable  owing  to  power  failures,  and  data  were  recorded  on  analog  tape  using  a 
sequence  of  seven  different  frequencies  ranging  from  15  to  45  kHz  in  5  kHz  steps.  A 
1  ms  pulse  length  was  used,  and  the  time  between  pulses  of  different  frequencies  was 
100  ms.  Each  sequence  of  pulses  was  transmitted  once  every  1.5  s,  for  a  total  averaging 
time  of  150  s. 

To  obtain  the  average  scattered  intensity,  a  thresholding  technique  was  devised  to 
align  the  pings  in  time.  Thresholding  for  the  forward  scattered  signal  was  performed  on 
the  combined  direct  and  single  bottom  bounce  paths  (time  difference  less  than  0.1  ms 
between  paths).  This  removed  travel  time  variations  due  to  motion  of  the  transmitter  in 
the  tidal  currents.  The  strong  air-sea  interface  return  was  used  for  thresholding  in  the 
vertical  incidence  measurements.  Thus  the  scattering  strength  profile  could  be  deter¬ 
mined  as  a  function  of  depth  below  the  instantaneous  sea  surface,  rather  than  the  mean 
sea  surface.  When  the  bubble  density  and  receiver  gain  were  high  (thus  saturating  the 
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air-sea  interface  return),  it  was  difficult  to  implement  this  thresholding  technique  on  the 
vertical  incidence  signal.  The  difficulty  arose  when  the  bubble  signal  level  became  com¬ 
parable  to  the  interface  signal  level,  and  the  bubble  signal,  rather  than  the  interface  sig¬ 
nal,  triggered  the  threshold.  In  such  cases  a  large  proportion  of  individual  pings  had  to 
be  carefully  examined  to  determine  the  proper  temporal  shift. 
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IV.  EXPERIMENTAL  RESULTS 

The  vertical  incidence  backscatter  results  will  be  presented  first,  followed  by  the 
forward  scattering  loss  measurements.  Table  1  summarizes  all  the  data  and  the 
corresponding  environmental  conditions  in  chronological  order. 


Table  1.  Acoustic  data  parameters  and  environmental  conditions. 


Run 

No. 

Frequency 

<\Hi) 

<7/ 

(dB) 

Oo.ffi 
(dB  re  m'1) 

L0.L, 

<m) 

?  B 

(m) 

SB  L 
(dB) 

Wind 

Speed 

(ms1) 

Wind 

Direcuon 

rms  Wave 
Height 
(m) 

43 

30 

.19  2 

-27.0 

0.61 

13 

W 

0.34 

44 

30 

-294 

-344,-39  1 

6.3,0.61 

4.27 

... 

13 

W 

0.40 

47 

30 

■236'’ 

... 

... 

23.1 

13 

W 

48 

30 

-24.6 

-22.5,-38  0,-40.5 

0.61.6.3.0.61 

2.05,6.2 

... 

13 

w 

050 

57 

25 

-21.0 

-25.0.-34.6 

2.49,1.26 

5.1 

13.5 

w 

066 

60 

25 

-25.2- 

... 

... 

15.0 

13.5 

w 

62 

25 

-24.8 

23.9 

0.82 

... 

13.5 

w 

0  51 

67 

20 

-34  3 

-30.1 

045 

10.8 

w 

044 

68 

20 

•26  0“ 

... 

... 

... 

11.5 

108 

w 

72 

20 

.37.9 

-38  1 

1  05 

10.8 

sw 

0.38 

75 

30 

•37.5 

-37  6.-40.5 

2  67,0.39 

1.32 

13 

sw 

0.36 

78 

40 

•32  4 

■32.1. -44  8 

1  00.0.38 

267 

... 

13 

sw 

0.47 

83 

20 

■52.9 

-49.2 

036 

7,6 

SE 

0.10 

84 

30 

-42  1 

-403 

066 

7.6 

SF. 

0.10 

88 

40 

-40  1 

-40  1.-44  7 

1  79.0.42 

1.48 

7.6 

SE 

0.10 

89 

sn 

-38.8 

-340 

0.57 

7.6 

SF. 

0 10 

91 

50 

05 

4.5 

SF. 

0.02 

134 

so 

■25  3 

■29. 5,-41  4 

3  10,0  57 

5  82 

120 

SE 

017 

13? 

20 

•30  1 

-29  0 

0  79 

120 

SE 

017 

136 

40 

-252 

-26  5. -32  9 

1  91,0  52 

2  35 

120 

SF. 

0.17 

137 

30 

•31  1 

-29  4 

067 

120 

SF. 

0.17 

138 

25 

-23  5 

•19  1 

0  46 

120 

SE 

0.17 

139 

50 

-24  0 

■24  5. -406 

1  21,061 

3  16 

120 

SE 

0.17 

32  1 

20 

■37  1 

■33  6 

0  44 

124 

SF. 

0.18 

3M 

25 

26  5 

•24  6.-44  0 

0  69.0  33 

2  19 

12  4 

SE 

0  18 

321 

30 

•36  1 

-34  1 

0  63 

124 

SF. 

0.18 

321 

35 

-32  1 

•26  3 

026 

124 

SE 

0.18 

32 

40 

•310 

-30  9 

0  97 

12  4 

SE 

0.18 

32-8 

20 

-38  4 

35  2 

051 

12  4 

SE 

0.18 

32  8 

25 

-27  7 

-22  4 

0  30 

12  4 

SE 

0.18 

32  8 

30 

-35  8 

•31  4 

0  37 

124 

SE 

0.18 

32-8 

35 

-31.2 

-27  3 

042 

12  4 

SK 

0.18 

32-8 

40 

-31  2 

-27  9 

0  47 

12  4 

SE 

0  18 

32-8 

45 

-27  0 

-24  8 

0  61 

124 

SE 

0.18 

32  3 

15 

-29  8" 

5.5 

12  4 

SE 

32  3 

20 

-29  O'* 

7.3 

12  4 

SE 

32-3 

25 

-28  8" 

8.2 

12  4 

SE 

32-3 

30 

•28  0" 

105 

124 

SF. 

32  3 

35 

•26  8" 

14  4 

12  4 

SE 

32  3 

40 

-267- 

12  7 

12  4 

SF 

32  10 

20 

•  30  9* 

124 

SE 

32  10 

30 

■30  2‘ 

124 

SE 

32-10 

40 

28  71, 

124 

SF. 

31  10 

20 

-44  7 

-40  0 

0  35 

9.7 

E 

0  10 

31  10 

2? 

•43  9 

-397 

0  38 

92 

E 

0 10 

3110 

30 

-41  4 

•  396.  470 

081.0  26 

1  42 

92 

E 

0 10 

31  10 

35 

-38  9 

3?  7 

0  48 

92 

E 

0 10 

31  10 

40 

-39  8 

•  38  4 

0  74 

9  2 

F. 

0 10 

31  10 

45 

•34  5 

-35  6 

1  31 

92 

F. 

0  10 

31-12 

20 

■  30  9" 

4  7 

92 

E 

71  12 

30 

-30  2* 

6  3 

92 

E 

3112 

40 

28  7" 

99 

92 

F. 

a  Inferred  O/  using  surface  bubble  loss 
^Inferred  O/  using  low  granng  angle  backscatter 
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A.  Vertical  Incidence  Backscatter 

The  average  volume  scattering  strengths  at  50  kHz  and  20  kHz  and  two  wind  speed 
conditions  (2  and  12  m/s)  are  shown  in  Figure  3  as  a  function  of  travel  time.  This  travel 
time  is  equivalent  to  depth  from  the  instantaneous  sea  surface  when  referenced  to  the 
air-sea  interface  return.  Note  that  when  the  conditions  are  nearly  calm  no  discernible 
bubble  return  is  observed,  and  the  leading  edge  of  the  air-sea  interface  return  is  sharp, 
rising  40  dB  in  about  1/3  ms.  This  clearly  indicates  that  bubbles  produced  by  mechan¬ 
isms  other  than  wave  breaking  (MacIntyre,  1986)  are  negligible.  This  simplifies  the  data 
analysis  and  allows  us  to  focus  solely  on  bubbles  created  by  wave  breaking  and  their 
acoustical  effects.  At  a  relatively  high  wind  speed  of  12  m/s  the  bubble  returns  are 
plainly  visible  above  the  background  signal  level  at  both  20  kHz  and  50  kHz.  The 
expected  exponential  rise  in  scattering  level  can  be  seen  at  20  kHz  from  the  straight  line 
visual  fit  to  the  data.  In  the  50  kHz  data  there  is  a  significant  change  in  slope  of  this  line; 
e.g.,  in  the  top  line  the  fitted  exponential  depth  constant  decreases  from  3.10  to  0.57  m. 
Thus  at  50  kHz  the  behavior  of  the  differential  backscattering  cross  section  per  unit 
volume  o5  due  to  bubbles  as  a  function  of  depth  is  described  as  follows; 

_  _  °oe~2'Lo  ,  0  <2<zb 

~  -zlL  0) 

o{e  1  ,  z b<z<°°  , 

where  z  is  positive  downward  from  the  sea  surface,  and  zg  represents  a  break  in  the 
exponential  decay  profile.  The  other  variables  in  Eq.  (1)  are  discussed  further  in 
Appendix  A.  Most  often  the  data  evidenced  an  exponential  profile  (see  Table  1;  entries 
with  a  single  depth  constant,  L0),  but  two-segment  profiles  like  those  shown  in  Figure  3 
were  relatively  common,  and  in  one  case  a  three-segment  profile  was  observed.  The 
deviation  from  an  exponential  profile  tended  to  occur  more  often  at  higher  wind  speeds 
and  frequencies.  We  note  that  the  exponential  profile  is  consistent  with  previously 
reported  trends  in  bubble  population  data,  e.g.,  Wu  (1981)  and  Kerman  (1986). 

A  convenient  and  informative  method  for  condensing  the  data  for  each  run  to  a  sin¬ 
gle  number  is  to  compute  the  integrated  volume  scattering  strength.  This  dimensionless 
quantity  is  found  by  integrating  under  the  bubble  return  portion  of  the  average  profile 
shown  in  Figure  3.  An  example  of  this  calculation  is  given  in  Appendix  A.  Figures  4 
and  5  are  examples  of  the  integrated  volume  scattering  strength  covering  the  frequency 
range  from  20-50  kHz  and  nominal  wind  speed  range  from  7.6-12.4  m/s  (as  measured 
from  Smith  Island).  The  expected  increase  in  scattering  level  with  wind  speed  is  clearly 
evident.  Data  taken  at  wind  speeds  less  than  3  m/s  showed  no  measurable  bubbles. 
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Figure  3.  T/  'ume  scattering  strength  averaged  over  50  pings  (200  s)  versus  time  into 
i,.j  ping  cycle.  The  individual  pings  are  aligned  in  time  with  respect  to  the 
instantaneous  sea  surface  by  thresholding  on  the  surface  returns.  Three  data 
runs  at  50  kHz  are  shown  as  light  solid  lines,  and  one  data  run  at  20  kHz  is 
shown  as  a  dashed  line.  The  straight  lines  represent  visual  fits  to  the  data. 
Bubbles  are  clearly  visible  for  the  two  50  kHz  runs  and  one  20  kHz  run  taken 
at  12  mis  ( labeled  Bubble  Returns);  they  are  not  visible  on  the  50  kHz  run 
taken  under  nearly  calm  conditions  (2  m/s).  The  nearly  calm  data  have  been 
adjusted  down  15  dB. 

The  frequency  dependence  of  the  integrated  volume  scattering  strength  is  also  of 
interest.  However,  we  expect  C/  to  be  a  strong  function  of  wind  speed,  and  therefore 
small  variations  in  U  during  the  course  of  an  experimental  run  can  mask  any  trend  with 
frequency.  Obtaining  the  data  in  Figure  4  during  the  early  pan  of  the  experiment  took 
approximately  an  hour  to  cover  the  frequency  range.  Because  those  in  Figure  5  were 
obtained  during  the  latter  pan  of  the  experiment  (in  the  multifrequenc’  mode),  the  same 
frequency  range  was  covered  in  about  5  minutes.  Thus  we  focus  c  n  Figure  5,  which 
shows  two  runs  taken  at  12.4  m/s  wind  speed  (separated  in  time  by  50  minutes)  and  one 
run  taken  at  9.2  m/s.  Although  there  is  only  one  value  available  at  45  kHz,  the  two  runs 
taken  at  12.4  m/s  are  extremely  consistent.  The  solid  lines  in  the  figure  represent  a 
power-law  fit  to  the  functional  form 

o/  =  afm  (2) 
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using  the  method  of  least  squares.  For  the  9.2  m/s  data  the  estimate  for  the  parameter  m 
is  2.60  (±1.5);  for  the  12.4  m/s  data  it  is  2.8  (±1.9).  We  report  the  optimal  values  for  the 
parameter  estimates,  with  the  bounds  in  parentheses  representing  the  95%  confidence 
interval  based  on  assumptions  regarding  the  underlying  distribution  of  the  data  (Bendat 
and  Piersol,  1986).  The  25  kHz  data  (Table  1,  runs  32-8  and  32-1)  were  not  included  in 
the  least  squares  analysis.  We  are  not  convinced  that  the  unusually  high  values  for  a/  at 
25  kHz,  measured  during  these  runs  as  well  as  some  earlier  runs,  are  accurate;  experi¬ 
mental  records  show  an  intermittent  interference  (of  unknown  source)  centered  near  the 
25  kHz  band. 


Figure  4. 


z  I 

LkJ 

cc 

►— 

(/)  O 

o  ° 

Z  i 

cr 


< 

o 

00 


<  ' 
ct:  ’ 
o 


z  ° 

o 
—  lO 


- 1 

I 

1 

- 1 - 

□  =  12 
0=76 

m/s 

m/s 

□ 

o 

□ 

- 

□ 

□ 

- 

- 

o 

0 

O 

- 

o 

_ t  —  ■- 

— -1 - 

-  -I.—, 

_ I _ 

10.0  20  0  30  0  40  0  50  0  60  0 

FREQUENCY  (kHz) 


The  integrated  volume  scattering  strength  attributable  to  bubbles  versus 
transmit  frequency,  at  two  wind  speeds. 


Figures.  The  integrated  volume  scattering  strength  attributable  to  bubbles  versus 
transmit  frequency.  The  upper  set  of  data  was  taken  at  a  wind  speed  of  12.4 
m/s,  and  the  lower  set  was  taken  at  9.2  m/s.  The  solid  lines  represent  a  least 
squares  fit  of  Eq.  (2)  to  the  data. 
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Limited  averaging  time  for  the  50-ping  runs  can  also  lead  to  variations  in  the  aver¬ 
age  O]  values,  particularly  since  the  spatial  sampling  size  was  relatively  small  (surface 
patch  size  of  about  1  m  diameter).  Figure  6  is  an  example  of  the  variation  in  the 
integrated  scattering  strength  as  a  function  of  time,  taken  during  the  early,  single  fre¬ 
quency  part  of  the  experiment.  The  interesting  characteristic  of  the  time  variation  of  C/ 
is  that  at  the  higher  wind  speeds  there  appear  to  be  longer  time  scales  of  variation.  Thus 
a  longer  time  averaging  period  is  necessary  if  the  variance  associated  with  the  final  mean 
value  is  to  be  reduced.  The  time  for  each  run  made  at  a  single  frequency  cannot  be 
extended,  however,  without  risking  significant  changes  in  the  environmental  conditions. 
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B.  Surface  Forward  Scattering 

Surface  forward  scattering  measurements  at  low  angles  were  interleaved  between 
the  vertical  incidence  measurements.  Figure  7  shows  two  examples  of  the  forward  scat¬ 
tered  signal,  one  for  low  wind  speed  and  one  for  a  relatively  high  wind  speed.  The  com¬ 
bined  direct  and  bottom  bounce  path  signals  are  evident  as  the  early,  near-replica  arrival, 
which  is  followed  by  the  surface  bounce  and  other  multipath  arrivals.  The  first  surface 
bounce  arrivals  rapidly  increase  to  a  peak  over  a  one-pulse-length  period  and  then  gradu¬ 
ally  decay  beyond.  With  the  low  wind  speed,  another  peak  occurs  at  12  ms  following  the 
direct  path  arrival,  which  is  due  to  the  multipaths  involving  two  surface  bounces.  Note 
that  the  curve  for  the  higher  wind  speed  is  about  10  dB  lower  in  level  and  displays  a 
much  slower  decay  in  level  beyond  the  first  peak;  also,  no  second  peak  is  visible. 


00  5.0  10.0  15  0  20.0  25.0  30.0  35.0  40.0 


TIME  (ms) 


Figure  7.  Normalized  forward  scattering  level  averaged  over  50  pings  (200  s)  versus 
time  into  the  ping  cycle.  The  normalization  is  such  that  a  single  direct  path 
should  appear  as  a  transmit  pulse  replica  with  a  level  of  0  dB  if  spherical 
spreading  plus  chemical  absorption  describes  the  path  loss.  The  upper  set  of 
data  was  taken  at  50  kHz  with  a  wind  speed  of  45  m/s,  and  the  lower  set  at 
20  kHz  and  wind  speed  of  10.8  m/s. 

The  forward  scattering  loss  due  to  a  single  surface  bounce  was  computed  by 
integrating  the  scattered  intensity  from  the  leading  edge  of  the  single-surface  bounce 
arrivals  to  the  leading  edge  of  the  double-surface  bounce  arrivals.  From  the  decibel 
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quivalent  of  this  integral,  1.5  dB  is  subtracted  to  account  for  the  inclusion  of  arrivals 
involving  bottom  bounces.  The  1.5  dB  value  is  determined  from  the  integrated  intensity 
measured  during  calm  conditions.  In  the  absence  of  bubble  attenuation  the  theoretical 
value  for  this  loss  is  0  dB;  i.e.,  no  loss  is  associated  with  scattering  from  the  air-sea  inter¬ 
face,  and  scattering  is  treated  as  a  0  dB  reflection  loss  (Thorsos,  1984a).  Any  measured 
loss  is  then  attributed  to  bubble  attenuation  effects  and  is  referred  to  as  surface  bubble 
loss  (SBL).  For  low  wind  speed  the  SBL  was  0.5  dB  (Table  1,  run  91).  For  high  wind 
speed  the  SBL  was  1 1.5  dB  (Table  1,  run  68). 

The  integrated  volume  scattering  strength  computed  from  the  vertical  incidence 
backscatter  measurements  can  be  compared  with  this  same  quantity  inferred  from  the 
surface  forward  scattering  measurements.  The  inferred  O/  is  calculated  from  the  SBL 
estimates  via  the  transfer  function 

cinft  Sp 

O/  [dimensionless]  =  —SBL  [dB]  ,  (3) 

109.2  6 


where 


5  =  2.55  x  10  3  f1/3  (fit  to  data  in  Devin,  1959) 
(total  damping  coefficient) , 

8*  =  0.0136 

(radiation  damping  coefficient), 


0  =  6.9° 

(grazing  angle  for  single  surface  bounce). 

This  calculation  is  reviewed  in  Appendix  B. 

Figures  8  and  9  compare  O/  derived  from  direct  vertical  incidence  measurements 
(VI)  with  that  inferred  from  surface  forward  scattering  loss  measurements  (FL)  taken  a 
few  minutes  apart.  In  Figure  9  we  have  also  plotted  estimates  of  C/  derived  from  low- 
angle  backscatter  measurements  for  0  =  30°  grazing  angle,  where  backscattering  strength 
is  convened  to  O/  following  the  method  described  by  McDaniel  (1988).  The  VI  meas¬ 
urements  in  Figure  9  are  an  average  of  those  shown  in  Figure  5,  which  were  taken  before 
and  after  the  FL  measurements.  The  solid  lines  represent  a  least  squares  fit  of  the  data  to 
Eq.  (2).  In  this  case  m  =  0.7*  for  the  FL  measurements  taken  at  U  =  9.2  m/s  (Figure  8), 
and  m  -  0.9  (±0.5)  for  the  FL  measurements  taken  at  U  =  12.4  m/s  (Figure  9). 


*The  three  FL  data  points  in  Figure  8  are  not  enough  to  compute  a  meaningful 
confidence  interval  about  the  optimum  estimate. 
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o 


Figure  8.  A  comparison  of  the  integrated  volume  scattering  strength  obtained  directly 
from  vertical  incidence  backscattering  measurements  (VI)  with  those  inferred 
from  surface  forward  scattering  loss  measurements  (FL).  The  data  were 
taken  at  a  wind  speed  of  92  mis.  The  solid  lines  represent  a  least  squares  fit 
of  Eq.  (2 )  to  the  data. 


O 


Figure  9.  A  comparison  of  the  integrated  volume  scattering  strength  obtained  directly 
from  vertical  incidence  backscattering  measurements  (VI)  with  those  inferred 
from  surface  forward  scattering  loss  measurements  (FL),  and  low-angle 
backscattering  measurements.  The  data  were  taken  at  a  wind  speed  of 
12.4  mis.  The  solid  lines  represent  a  least  squares  fit  of  Eq.  (2)  to  the  VI  and 
FL  data. 
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V.  DISCUSSION  AND  SUMMARY 

There  is  a  significant  offset  between  the  VI  and  FL  estimates  of  G/,  but  the  FL  and 
low  angle  backscatter-derived  estimates  are  within  experimental  error  of  each  other.  The 
reason  for  this  difference  may  lie  in  the  existence  of  a  thin  (-  10  cm)  layer  of  high  bubble 
concentration  next  to  the  air-sea  interface.  Such  a  near-surface  layer  has  been  observed 
and  is  termed  the  "bubble  generation  layer"  (Baldy,  1988).  The  finite  pulse  lengths  used 
in  the  VI  measurements  could  not  detect  bubbles  within  such  a  thin  layer  and  resolve 
them  from  the  air-sea  interface  return — despite  our  attempts  to  correct  for  this,  as  dis¬ 
cussed  in  Appendix  A.  On  the  other  hand,  both  the  FL  and  low-angle  backscatter  meas¬ 
urement  techniques  inherently  sample  this  layer.  An  additional  observation  lending 
further  support  to  the  existence  of  such  a  layer  is  that  the  VI  data  slowly  converge  toward 
the  FL  data  at  the  higher  frequencies  corresponding  to  smaller  resonant  bubble  radii. 
These  smaller  bubbles  have  lower  buoyancy  rise  rates,  and  thus  have  longer  residence 
times  below  the  sea  surface  and  within  a  zone  that  can  be  resolved  by  the  VI  measure¬ 
ments.  We  also  note  that  the  smaller  difference  between  the  VI  data  and  the  FL  data  at 
the  higher  of  the  two  wind  speeds  is  suggestive  of  the  downward  displacement  of  a 
greater  proportion  of  the  bubbles  as  turbulent  mixing  increases  with  wind  speed. 

Returning  to  Figures  8  and  9,  we  note  the  difference  in  the  frequency  power  laws 
for  CT/  derived  from  VI  and  FL  data;  this  difference  is  approximately  1.9  in  the  exponent 
m  in  Eq.  (2)  for  both  figures.  An  equivalent  frequency  power-law  relation  for  NR ,  the 
depth  integrated  bubble  density  at  resonance,  can  be  derived  from  Appendix  B.  The  VI 
data  show  a  frequency  dependence  for  N R  of  approximately  /6;  for  the  FL  data  this 
dependence  is  approximately  y .  The  frequency  power  law  derived  from  FL  data  is  con¬ 
sistent  with  the  data  of  Farmer  and  Lemon  (1984)  taken  at  14.5  kHz  and  25  kHz.  Furth¬ 
ermore,  the  frequency  power  law  derived  from  VI  data  suggests  that  NR  falls  off  more 
rapidly  with  size  than  would  be  indicated  by  the  FL  data.  This  is  also  consistent  with  the 
VI  approach  not  being  sufficiently  sensitive  to  the  bubble  population  in  the  near-surface 
layer. 

Finally,  there  is  a  significant  difference  in  the  wind  speed  dependence  for  the  fitted 
curves  shown  in  Figures  8  and  9.  For  example,  for  the  wind  speed  change  from  9.2  to 
12.4  m/s,  both  the  VI  and  FL  fitted  curves  show  a  nearly  constant  displacement  across  the 
frequency  range  20-40  kHz:  about  6  dB  for  the  VI  data  and  about  2  dB  for  the  FL  data. 
In  order  to  infer  a  wind  speed  power-law  dependence  while  accounting  for  the  frequency 
dependence,  and  thereby  utilize  all  the  data  in  Table  1 ,  we  proceed  with  a  multivariable 
power-law  fit  to  the  functional  form 

o,  =  &fmUn  (4) 

using  the  method  of  least  squares.  A  new  arbitrary  constant  p  is  used  here  to  distinguish 
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it  from  the  a  in  Eq.  (2).  The  optimal  estimate  for  the  parameter  m  is  2.4  (±0.6)  for  the  VI 
data  and  0.6  (±0.4)  for  the  FL  data.  In  both  cases  the  estimates  lie  within  the  confidence 
intervals  computed  using  only  the  multifrequency  data  of  Table  1  (from  run  32-1  down), 
showing  the  consistency  in  the  entire  data  set.  The  optimal  estimate  for  the  parameter  n 
is  6.1  (±1.6)  for  the  VI  data,  which  indicates  O/  «  £/61,  or  equivalently,  Nr  «  U6A,  and 
is  consistent  with  data  from  Farmer  and  Lemon.  On  the  other  hand,  for  the  FL  data  the 
optimal  estimate  for  n  is  2.2  (±1.6),  indicating  O/ «  U22,  or  equivalently,  Nr  <*U22. 
However,  other  reports  (e.g.,  Thorsos,  1984b)  have  shown  surface  forward  scattering 
data  similar  to  ours  that  are  not  uniquely  determined  by  wind  speed.  Furthermore,  com¬ 
parisons  and  interpretations  of  the  wind  speed  and  frequency  power  laws  reported  here 
are  tenuous  because  the  fetch-limited  conditions  off  Whidbey  Island  differ  from  those  in 
the  open  ocean,  e.g.,  in  the  length  scales  of  sea  surface  waves  and  characteristics  of  bub¬ 
ble  plumes. 

In  summary,  it  is  clear  from  the  vertical  incidence  measurements  that  bubbles  are 
acoustically  observable  at  wind  speeds  as  low  as  3  m/s  and  that  for  wind  speeds  greater 
than  about  5-6  m/s  the  surface  bubble  loss  (SBL)  can  become  quite  large  (>10dB). 
Direct  estimates  of  the  depth  integrated  volume  scattering  strength  O/  were  obtained 
from  vertical  incidence  measurements  (VI  data),  and  inferred  estimates  of  O/  were 
obtained  from  the  SBL  estimates  associated  with  the  low  angle  surface  forward  scatter¬ 
ing  loss  measurements  (FL  data).  In  most  cases  the  VI  and  FL  measurements  were 
separated  in  time  by  no  more  than  a  few  minutes.  The  principal  results  are 

(1)  The  depth  dependence  of  the  near-surface  volume  scattering  strength  fol¬ 
lowed  the  expected  exponential  behavior  for  the  majority  of  cases,  but 
significant  departures  occurred  in  some  cases  where  a  two-part  exponential 
behavior  is  observed. 

(2)  The  frequency  dependence  of  a /  derived  from  VI  data  was  a/  «  f2A\  the  fre¬ 
quency  dependence  of  Cj  inferred  from  FL  data  was  O/  «  /°-6. 

(3)  The  wind  speed  dependence  of  0/  derived  from  VI  data  was  0/  «  U6A ;  the 

wind  speed  dependence  of  C/  inferred  from  FL  data  was  0/  U22. 

(4)  The  differences  in  the  VI  and  FL  data  noted  in  (2)  and  (3)  may  in  part  be  due 
to  the  VI  data  being  less  sensitive  than  the  FL  data  to  the  near-surface  bubble 
layer. 
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APPENDIX  A 

This  appendix  discusses  the  variables  in  Eq.  (1)  and  reviews  the  calculations 
involved  in  obtaining  an  estimate  of  Gj.  Figure  A1  shows  a  profile  of  the  reverberation 
level  (in  arbitrary  decibel  units)  as  a  function  of  time.  The  profile  rises  exponentially 
above  the  noise  level  between  the  points  marked  by  2  and  1,  owing  to  scattering  from 
bubbles.  We  first  determine  an  exponential  depth  constant  L0  that  governs  the  profile 
between  points  1  and  2,  such  that 

h=  I\e~2'L°  .  (Al) 

For  example,  for  the  levels  shown  in  Figure  Al  we  find  La=  0.303  m. 


Figure  Al.  Profile  of  the  reverberation  level  (in  arbitrary  decibel  units)  averaged  over 
50  pings  versus  time  into  the  ping  cycle.  The  behavior  between  points  1  and 
2  is  modeled  as  an  exponential  decay. 

Next  I0,  the  maximum  intensity  in  the  reverberation  profile,  is  determined.  By 
hypothesis,  I0  is  due  to  backscattering  from  bubbles  located  immediately  below  the 
surface  but  not  from  the  air-sea  interface  itself.  That  is,  we  remove  the  effects  of  sea 
surface  backscatter  in  the  midst  of  a  reverberation  producing  volume.  However,  since 
the  reverberation  profile  is  dominated  by  the  air-sea  interface  beyond  point  1,  we  infer  a 
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level  for  I0  by  extending  the  exponential  profile  a  distance  ct/4  =  0.364  m  beyond  point 
1.  The  distance  ct/4  is  used  because  it  centers  the  equivalent  scattering  volume 
(producing  10)  at  the  air-sea  interface.  For  this  case  we  find  1Q  =  3.6  dB.  By  including 
system  gains  and  receiving  sensitivity,  this  is  equivalent  to  a  reverberation  level  RL0  = 
120dB  re  | xPa. 

The  reverberation  profile,  which  has  now  been  extended  to  the  air-sea  interface,  is 
converted  to  volume  scattering  strength  Sy0  in  a  straightforward  manner  following  Urick 
(1983).  For  example, 

SVo  =  RL0  -  SL  +  40  log  R-  10  log  V  .  (A2) 

The  relevant  parameters  for  this  example  are  SL  =  204.8  dB  re  \xPa  and  R  =  29.8  m,  with 
the  scattering  volume  V  =  (cx/2)R2\\f,  where  \y  is  the  equivalent  solid  angle 
corresponding  to  the  combined  transmitting  and  receiving  line  arrays.  We  account  for 
sound  absorption  in  seawater  and  arrive  at  an  initial  estimate  of  the  volume  scattering 
strength  for  the  near-surface  layer  due  to  bubbles.  We  refer  to  this  as  S  Vo.  which  is 
equal  to  -23.4  dB  re  1  m-1 ,  or  equivalently  a '0  -  0.00457. 

A  final  estimate  of  a0  and  (T/  must  also  include  the  effect  of  absorption  due  to 
bubbles.  The  two-way  transmission  loss  due  to  bubbles  is  expressed  as  (Medwin,  1977) 

TLb  =  8.686  /  ae(z)dz  ,  (A3) 

O 

where  oe  is  the  total  extinction  cross  section  per  unit  volume.  The  extinction  and 
differential  backscattering  cross  sections  per  unit  volume  due  to  bubbles  are  related  by 


where  8  and  bp  are  defined  below  Eq.  (3).  The  depth  dependence  of  a*  usually  assumes 

_  2  // . 

a  single  profile  equal  to  a0  (e  ° ).  In  some  cases  the  profile  may  consist  of  two  parts 
involving  Oj  and  L  \  as  shown  in  Eq.  (1).  We  use  the  single  profile  case  in  the  integral  in 
Eq.  (A3)  because  oQ  is  always  much  greater  than  Oj ,  giving 

TLb  =  8026  bo0L0  .  (A5) 

Finally,  an  equation  for  a  revised  estimate  of  Oa  that  includes  the  effect  of  bubble 
absorption  is 

10  log  <5a  =  5  Vo  +  8026  b  <5 o  La  .  (A6) 
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Equation  (A6)  has  only  one  unknown,  oot  which  we  readily  solve  numerically  using  the 
Newton-Raphson  method  (Hildebrand,  1976).  The  result  for  this  example  is  SVo  = 
-22.4  dB  re  1  m-1.  Finally,  <5/  is  computed  by  depth  integrating  a*;  e.g.,  in  the  single 
profile  case  we  have 

Oj  =  o0L0  .  (A7) 
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APPENDIX  B 

This  appendix  reviews  the  means  of  conversion  between  SBL  and  O/.  The 
backscattering  cross  section  per  unit  volume  due  to  entrained  bubbles  is  (Medwin,  1977) 

'X  A  ir  t 

Ar  (z)  N(Ar) / 8  ,  (Bl) 

4 

where  Ar  is  the  radius  of  resonant  bubbles,  N  (Ar)  is  the  number  of  resonant  bubbles  per 
unit  volume  per  unit  bubble  radius  in  microns,  and  5  is  defined  below  Eq.  (3).  Medwin 
also  shows  that  Ar  (in  microns)  at  depth  z  (in  meters)  is  related  to  frequency  (in 
kilohertz)  by 

Ar  *  3-5— 3  Vl  +  0.1  z  .  (B2) 

By  definition, 


o,  =  ~ir  i  (B3) 

O 

Because  Ar  is  slowly  varying  over  the  range  of  significant  bubble  concentration,  we  have 

oo 

a/*—15-  J  4  J  N(Ar)  dz  .  (B4) 

O 

The  integral  in  Eq.  (B4)  is  defined  as  Nr,  the  depth-integrated  bubble  density  at 
resonance.  SBL  is  defined  by  Thorsos  (1984a)  as 

oo 

SBL  (dB)  =  — ~ -  J  <5e{z)dz  .  (B5) 

sin6  ^ 


In  our  case  we  integrate  from  d  =0.  We  note  that  SBL  is  just  the  two-way  transmission 
loss  due  to  bubbles,  (Eq.  A3),  modified  for  arbitrary  angle  of  incidence.  From  Eqs.  (A4) 
and  (Bl),  and  carrying  out  the  integration  as  in  (B3),  we  find 

8.686  3.4:t2 


SBL  (dB)  = 


sinG  A3r(z)6r 


Hr, 


(B6) 


which  leads  to  the  transfer  function,  Eq.  (3)  in  Section  IV  of  the  text. 
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